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ABSTRACT 

The integrated colors of distant galaxies provide a means for interpreting the properties of their stellar con- 
tent. Here, we use rest-frame UV-to-optical colors to constrain the spectral-energy distributions and stellar 
populations of color-selected, B-dropout galaxies at z ~ 4 in the Great Observatories Origins Deep Survey. 
We combine the ACS data with ground-based near-infrared images, which extend the coverage of galaxies at 
Z ~ 4 to the rest-frame B-band. We observe a color-magnitude trend in the rest-frame m(UV)-B versus B 
diagram for the z ~ 4 galaxies that has a fairly well-defined "blue-envelope", and is strikingly similar to that 
of color-selected, {/-dropout galaxies at z ~ 3. We also find that although the co-moving luminosity density 
at rest-frame UV wavelengths (1600A) is roughly comparable at z ~ 3 and z ~ 4, the luminosity density at 
rest-frame optical wavelengths increases by about one-third from z ~ 4 to z ~ 3. Although the star-formation 
histories of individual galaxies may involve complex and stochastic events, the evolution in the global lumi- 
nosity density of the UV-bright galaxy population corresponds to an average star-formation history with a 
star-formation rate that is constant or increasing over these redshifts. This suggests that the evolution in the 
luminosity density corresponds to an increase in the stellar-mass density of > 33%. 

Subject headings: early universe — cosmology: observations — galaxies: evolution — galaxies: formation 
— galaxies: high-redshift — galaxies: photometry 



1. INTRODUCTION 

Current investigations of high-redshift (z > 2) galaxies 
have been focusing on the properties of these objects as 
a global population. Many surveys identify these galax- 
ies by their strong emission at rest-frame UV wavelengths 
(observed-frame optical) and spectral breaks at Lyman a and 
the Lyman limit (so-called Lyman-break galaxies [LBGs]; 
e.g., Giavalisco 2002). These galaxies are generally domi- 
nated by the light from OB stars, and ha ve properties that ar e 
similar to local starburst galaxies (e.g., Sha plev et alJl2003l) . 
Near-infrared (NIR) photometry of galaxies at z > 2 extends 
the observations to rest-frame optical wavelengths, probing 
the light from A- and later-type stars. Several studies have 
used NIR observations to constrain the properties of the stellar 
populations of z ~ 2-3 galaxies (e.g., Sawicki & Yee 1998; 
Papovich, Dickinson, & Ferguson 2001; Shapley et al. 2001; 
Labbe et al. 2003; Franx et al. 2003), and to estimate the evo- 
lution of the global stellar-mass density for < z < 3 (e.g., 
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Dick inson et all2003l) . 

At present, some of the constraints on the parameters of 
the galaxies' stellar-population models are uncertain by more 
than an order of magnitude (Papovich et al. 2001; Shapley 
et al. 2001). Even so, the stellar-population ages and star- 
formation histories of the models have broad implications for 
galaxy evolution at higher redshifts (see Ferguson, Dickin- 
son, & Papovich 2002). The galaxies' spectral-energy dis- 
tributions (SEDs) contain the integrated record of their past 
and current star formation. Thus, comparing galaxy SEDs at 
different redshifts allows us to improve the constraints on the 
star-formation histories of these galaxies. 

In this Letter, we study galaxies at z ~ 4 selected from 
deep imaging with the Advanced Camera for Surveys (ACS) 
onboard the Hubble Space Telescope {HST) and augmented 
with NIR observations from the ground (§ 2), and we com- 
pare these to similar rest-frame colors for color-selected 
galaxies at z ~ 3 from the Hubble Deep Fields, North and 
South (HDF-N and -S; § 3). We then discuss the SEDs 
of the luminosity density generated by these galaxies, and 
we consider the implications on the galaxies' star-formation 
histories (§ 4). Throughout this Letter, we use a flat cos- 
mology with il m = 0.3, Qa = 0.7, a Hubble constant of 
70 km s -1 Mpc -1 , and we use AB magnitudes, wab = —48.6 — 
2.51og(/ 1 ,/l erg s" 1 cm" 2 Hz" 1 ). 

2. THE OBSERVATIONS AND GALAXY SAMPLE 

The HST data used in this letter stem from the first three 
epochs of ACS imaging of the Chandra Deep Field South 
(CDF-S) as part of the Great Observatories Origins Deep Sur- 
vey (GOODS) program. These data provide imaging in the 
B435, V606> '775 > an d Z85o-bands covering a field of view of 
10' x 16.5'. The dataset, its reduc tion procedures, a nd object 
cataloging are described in Giavalisco et al. (2003a). 

The ACS data only probe rest-frame UV wavelengths for 
galaxies at z > 1.2. To study the rest-frame optical light 
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FIG. 1. — Rest-frame UV-optical color-magnitude diagrams for color-selected high-redshift galaxies. Left: Panel shows the i^;—K s versus K s color- 
magnitude diagram for the B-dropout galaxies selected from the GOODS CDF-S data (squares; triangles denote lcr limits). These filters correspond to rest- 
frame 1600- and 4400-A colors. The dashed line illustrates the approximate zgso-band detection limit (for (775 -Z850 ~ 0)- Right: Panel shows the V^—H 
versus H color-magnitude diagram for f/-dropout galaxies selected from the WFPC2 data for the HDF-N (filled circles) and HDF-S (open stars). The abscissa 
of each panel has been adjusted to show an approximately equal range of rest-frame B-band absolute magnitudes. The dotted line shows the magnitude limit 
(m* + 1 mag) used to define samples that are used to derive the luminosity densities in § 3. The solid line in each panel indicates the fiducial 'blue-envelope' . 



from higher redshift galaxies we combined the ACS data with 
VLT/ISAAC imagin g in the JHK S bands obtained as part of 
GOODS (Giavalisco et al. 2003a). Currently, we use eight in- 
dividual ISAAC tiles covering a total area of 50 arcmin . 

We have analyzed the ISAAC JHK S images to extract op- 
timal photometry matched to the ACS zgso-selected catalog. 
In brief, we use the zsso-band data to create two-dimensional 
templates for each object, convolve these to match the point- 
spread function (PSF) of the ground-based data, and scale 
and fit the convolved templates to the JHK s -band images to 
extract object fluxes (see C. Papovich et al. in preparation; 
Fernandez-Soto, Lanzetta, & Yahil 1999). The advantage of 
this method as applied to faint galaxies is that it reduces con- 
cerns about PSF- and aperture-matching effects on the rel- 
ative photometry and permits deblending of objects partially 
merged by the terrestrial seeing. 

We identified high-redshift galaxies whose SEDs reveal 
a spectral discontinuity at the observed Lyman limit, which 
arises from H I ab sorption systems along the line of sight 
(e.g.. lMadaulll995l) . We selected "fi-dropout" galaxies with, 
#435-V606 > 1 -1 + (Visnfi - Zsvi), - V fififi > 1.1, and V 606 - 
Z&50 < 1-6 (see Giavalisc o et al.ll2003bl) . Using these colors, 
we are able to reject most lower-redshift interlopers while de- 
tecting galaxies at z ~ 4 with bright rest-frame UV spectra 
(which is indicative of ongoing star formation; Leitherer et 
al. 1999). Indeed, using the ~ 1400 spectroscopic redshifts 
available for GOODS objects, we observe no contamination 
from stars or galaxies at z < 3. Our modeling indicates that 
these criteria identify galaxies with an expected redshift dis- 
tribution that peaks at z = 3.9 and tapers to higher and lower 
redshifts, with 50% completeness limits of 3.4 < z < 4.5. 

3. UV-OPTICAL COLORS OF HIGH-REDSHIFT GALAXIES 

In the left panel of figure we show the (775 - K s versus 
K s color-magnitude diagram for the Z?-dropout galaxies. For 
z ~ 4, these colors correspond approximately to rest-frame 



ra(1600 A) - m(4400 A). Although most of the fi-dropout 
galaxies have relatively blue rest-frame UV-to-optical col- 
ors, the brightest galaxies in K s have redder optical-infrared 
colors with a fairly well-defined "blue-envelope". Because 
this sample is selected in the zgso-band, there is no reason to 
expect that the lack of bright galaxies (K s < 24) with blue col- 
ors is due to some selection effect. At fainter magnitudes, part 
of this trend may be due to increasing photometric uncertain- 
ties. However, the trend among the B-dropouts is apparent 
even at bright magnitudes where the photometric errors are 
small, and exists for the {/-dropouts (see below), where the 
optical and IR data are much deeper than the GOODS data. 

A similar blue-envelope exists for the colors of galaxies at 
z ~ 3. The right panel of figure^shows the V(,q(,-H versus H- 
band color-magnitude diagram for {/-dropout galaxies from 
the WFPC2 and NICMOS data of the HDF-N (Papovich et 
al. 2001; M. Dickinson et al. in preparation) a nd from the 
WFPC 2 and VLT/ISAAC data of the HDF-S (Lab be et alJ 
120031) . and using the color-selection of Steidel et al. (1999). 
The {/-dropout criteria identify galaxies with a mean redshift, 
z ~ 2.6, and 50% completeness limits of 2.0 < z < 3.3. Al- 
though the galaxy samples for the HDF-N and -S are origi- 
nally detected in NIR data, our tests using an /-band-selected 
sample for the HDF-N yield essentially no change in the num- 
ber of selected {/-dropout galaxies nor in the number and 
luminosity densities derived below. The advantage of using 
these data is that they provide band-matched catalogs from 
optical-to-NIR wavelengths that are comparable to the data 
for the B-dropout galaxies. 

For the z ~ 3 galaxies, the V6O6 - H color corresponds to 
nearly identical rest-frame colors as 1775 - K s for the galax- 
ies at z ~ 4. The color-magnitude trends at z ~ 3 and 4 are 
strikingly similar. Thus, the LBGs at both z ~ 3 and z ~ 4 
seem to require that either (or likely some combination of) the 
mean stellar population age, metallicity, and/or dust opacity 
increase with optical luminosity. Furthermore, we also note 
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that the [/-dropout galaxies have slightly redder rest-frame 
m{ 1 600 A) - m (4400 A) colors as evidenced by the fact that the 
"blue-envelope" shifts by ~ 0.2 mag. This reddening seems 
to imply that the mean ages, metallicity, and/or dust opacity 
are actually increasing from z ~ 4 to 3. 

4. THE INTEGRATED SEDS OF HIGH— REDSHIFT GALAXIES 

To measure the number and luminosity densities, we re- 
quire the redshift distribution and effective volumes probed 
by the U— and B-dropout criteria. For the U -dropout galax- 
ies, we have used the effective volume derived by Steid el et alJ 
( 1999). For the B-dropout galaxies, we have simulated arti- 
ficial LBGs in the GOODS data with a distribution of col- 
ors and sizes that match the observed B-dropout properties 
(see Ferguson et al. 2003). We then measure photometry 
for the simulated galaxies and apply our color-selection cri- 
teria to derive the probability that a LBG with given mag- 
nitude and redshift is detected in the data, p(m,z). This 
technique accounts for incompleteness due to magnitude and 
surface-brightness effects as well as the color-selection pro- 
cess. We then compute the effective volumes, V e s(m) = 
J dz p(m,z)dV(z)/dz- Because these selection criteria are in- 
sensitive to heavily obscured galaxies and to galaxies with 
passively evolving, older stellar populations, our estimates 
of galaxy densities are strictly lower limits, derived from the 
UV-bright population alone. 

To calculate the galaxy densities we consider galaxies 
brighter than m* + 1 mag, where m*(3£) — 24.20 for the U— 
dropouts (Steidel et al. 1999; where 8% * (F 6 06 + 4m) A 
and m* has been adjusted slightly to account for differences 
in cosmology and the mean redshifts of the samples here 
and those of Steidel et al.), and m*(I) = 24.70 for the B- 
dropouts (Giavalisco et al. 2003b). We derive specific lu- 
minosity densities by integrating the flux densities in each 
bandpass to m* + 1 mag (where the B-dropout galaxies are 
well-detected in both the ACS and ISAAC data, see figure^: 
p v = JdL v (m)n(m)L v (m)/V e ft(m), where n(m) is the ob- 
served number density, L v (m) = 10- QA ^- 6+m Hndl(l +z)~\ 
and di is the luminosity distance. We plot the measured val- 
ues in figure [2] along with uncertainties from a bootstrap re- 
sampling. Note that we have made no correction for fainter 
galaxies. 

The area of the HDFs and the portion of the GOODS field 
with ISAAC imaging are fairly small, and cosmic variance 
may be a substantial source of uncertainty (see Somerville et 
al. 2003). To limit this effect on our results, we have normal- 
ized the luminosity densities of both the U— and B-dropout 
samples such that they match the rest-frame UV values of 
Steidel et al. (1999) and Giavalisco et al. (2003b), respec- 
tively, while preserving the measured "color" between the 
rest-frame UV and other wavelengths. These surveys aver- 
age over more area and sightlines than considered here, which 
significantly reduces the uncertainties from cosmic variance. 
Cosmic variance may still be a factor for the B-band lumi- 
nosity densities. The luminosity weighted mean rest-frame 
ra(1600A) -m(4400A) colors of the {/-dropout galaxies in 
the HDF-N and -S vary by ~ 0.15 mag, which provides an 
estimate on the uncertainty in the typical luminosity-density 
colors of LBGs in single HDF-sized fields. Future studies us- 
ing the entire GOODS dataset (with full {/-band and NIR ob- 
servations of both fields) will substantially increase the sam- 
ple size of both the U— and B-dropouts by roughly a factor of 
» 25 and 10, respectively, and should strengthen our results. 
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FIG. 2. — The co-moving luminosity density as a function of rest-frame 
wavelength. The data points show the values derived for the B-dropout galax- 
ies in the GOODS CDF-S data (filled squares, solid lines) and for the U— 
dropouts in the HDF-N and S (open stars, dashed lines). The luminosity den- 
sities have been integrated down to m* + 1 mag with no correction for galaxies 
fainter than this limit or due to dust extinction (see text). The shaded region 
spans the range of empirical starburst galaxy SEDs from Kinney et al] |19961) . 
from the bluest [NGC 1705, E(B-V) ~ 0.0] to reddest [0.6 < E(B-V) < 0.7] 
templates. The vertical dotted lines indicate Lyman a and the Balmer Break. 



The fact that the rest-frame 1600 A luminosity density is 
roughly unchanged (to within < 10%) implies that the star- 
formation rate (SFR; at least of massive, OB stars) is roughly 
constant (see, e.g., Madau, Pozzetti, & Dickinson 1998). 
However, we observe marginal evidence for a steepening of 
the UV continuum (~ 1500-2500 A), and strong evidence 
(w 98% confidence) that the average m(1600A) -m(4400A) 
color becomes redder for the {/-dropouts relative to those 
for the B-dropouts. In particular, the luminosity density at 
rest-frame 4400 A grows from 2.1 ± 0.3 x 10 26 erg s" 1 Hz" 1 
Mpc" 3 at z ~ 4 to 2.8 ± 0.3 x 10 26 erg s" 1 Hz" 1 Mpc" 3 at 
z ~ 3. Based on the observed color distribution of LBGs (to 
fainter magnitude limits than our samples here), our simula- 
tions suggest that this result is unlikely to be biased signifi- 
cantly by the color-selection of these galaxies. This change 
in the mean color suggests that the mean stellar mass-to-light 
ratio of these galaxies is also increasing. Therefore, the total 
stellar mass density of this population is presumably growing 
by more than the > 33% increase seen in the rest-frame B- 
band optical light. Dickinson et al. (2003) found evidence for 
a substantial build-up in the global stellar mass density from 
z ~ 3 to z ~ 1 . The present results extend that trend to still 
higher redshifts. 

The luminosity densities constrain the global, average star- 
formation histories of galaxies at z ~ 3-4. One possible sce- 
nario is a roughly constant SFR with an unchanging dust con- 
tent. This would produce approximately equal p„(1600A) at 
z ~ 3 and 4, and build up the global stellar mass with time, 
which would redden the average galaxy SED from z ~ 4 to 
z ~ 3 as the B-band light increases. We note that the Uni- 
verse ages by ~ 50% from z ~ 4 to 3, and may imply that 
galaxies formed stars at a near-constant rate from z 3> 6 to 3. 
An alternative scenario is a rising SFR from z ~ 4 to 3 with an 
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increasing dust content that is tuned to produce a roughly con- 
stant /^(l^OOA), and is consistent with the suggestion that the 
UV continuum is redder at z ~ 3. The fractional stellar-mass 
build-up in this scenario would be still greater than that with 
constant SFR and an unchanging dust content. 

Both of the above scenarios invoke a SFR that is either con- 
stant or growing in time. However, the stellar populations ob- 
served in z ~ 3 galaxies appear to be young (a few x 10 8 yr for 
simple, monotonic star-formation histories; Sawicki & Yee 
1998; Papovich et al. 2001; Shapley et al. 2001), which im- 
plies that these stars had not formed by z ~ 4. Galaxies at 
these redshifts probably have star-formation histories that are 
more complex (with multiple formation episodes) than those 
described by the simple models (see Ferguson et al. 2002). 
However, even though the star-formation histories for indi- 
vidual galaxies may be fairly stochastic, the global, average 
SFR at high redshift must be roughly constant in order to ex- 
plain the evolution in the UV and optical luminosity densities. 

Although the increase in the observed optical luminosity 
density from z ~ 4 to z ~ 3 appears robust, the uncertain- 
ties due to cosmic variance and other systematics are non- 
negligible. If there is, in fact, no evolution in the optical lumi- 
nosity density then this could imply that stars formed in the B- 
dropout galaxies are not present within the [/-dropout galax- 
ies, and that at some point the z ~ 4 galaxies cease forming 
stars and evolve beyond the [/-dropout selection criteria. This 
would have interesting implications for constraints on the 
sources of reionization and/or on the IMF at higher redshifts 
(see, e.g., Ferguson et al. 2002). However, this would likely 
be inc onsistent with the low number density of red, z > 2 ob- 
jects jFranx et alJl2003l) . Moreover, the majority of galaxy 
evolution models predict growth in the stellar-mass density 
by amounts consistent with our observations (although spe- 
cific predictions vary from model to model; see Dickinson et 
al. 2003, and references within), and this suggests that much 
of the stellar-mass assembly at high-redshifts does occur in 
these UV-bright galaxies. 

5. SUMMARY 

We have compared the rest-frame UV-optical colors of 
color-selected galaxies at z ~ 4 to those of similarly selected 



galaxies at z ~ 3. We find a great degree of similarity in the 
rest-frame m(UV)-B versus B diagram for galaxies from 
z ~ 3 - 4. However, although the rest-frame UV luminos- 
ity densities at z ~ 4 and z ~ 3 are comparable, there is ev- 
idence that the rest-frame fi-band luminosity density grows 
by w 33 ± 16%. Even though the star-formation histories of 
individual galaxies may involve complex and stochastic pro- 
cesses, the evolution of the luminosity density corresponds 
to a globally average SFR that is constant or increases with 
time. This implies that the average stellar-mass-to-light ratio 
of galaxies is also increasing over this redshift range and that 
the global stellar-mass density grows by more than > 33% 
over the ~ 1 Gyr interval that elapses between z ~ 4 and 3. 

By selecting in the rest-frame UV, we are likely to miss 
galaxies without ongoin g, relatively u nobscured star forma- 
tion. For example, Fra nx et al.l l|2£)03) have used deep NIR 
images of the HDF-S to identify red galaxies which may have 
evolved, massive stellar populations at z ~ 2, and may con- 
tribute significantly to the global stellar mass density. Such 
objects would be missing from the samples here, leading to an 
underestimate of the total density. Logically, the number and 
mass content of these red galaxies should grow as the Uni- 
verse ages. Therefore, the apparent increase in stellar mass 
density from z ~ 4 to 3, as traced here by the UV-bright pop- 
ulation, would only be strengthened if red, UV-faint objects 
were also considered. Ultimately, the full GOODS dataset 
will permit more direct tests of the star-formation histories of 
the distant galaxy population, as will future observations with 
the Space Infrared Telescope Facility and the James Webb 
Space Telescope. 
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